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BIRKHOFF'S THEOREM FOR PANSTOCHASTIC MATRICES 



Dean Alvis and Michael Kinyon 



1. Introduction. 

An n X n matrix witli nonnegative real entries is doubly stochastic if the sum of the 
entries along any of its rows or columns is equal to 1. A doubly stochastic matrix is 
panstochastic if the sum of the entries along any downward diagonal or upward diagonal, 
either broken or unbroken, is equal to 1. For example, the matrices 



1 

5 



11111 
11111 
11111 
11111 
11111 



1 

1 

10 

1 

1 OJ 



1 

60 



1 13 20 7 19 

22 9 16 3 10 

18 12 24 6 

14 21 8 15 2 

L 5 17 4 11 23 



are panstochastic. 

A linear combination is called convex if the coefficients are nonnegative and their sum 
is equal to 1. In [1], Birkhoff showed that every doubly stochastic matrix can be expressed 
as a convex combination of permutation matrices. Related results for integral matrices 
had been obtained earlier by Konig [8] and Egervary [5]. Birkhoff 's theorem has been 
generalized in various ways; for example, Schneider obtained the result for matrices with 
entries in lattice-ordered abelian groups [11]. 

Does the analogue of Birkhoff's theorem hold for panstochastic matrices? Our first 
main result is that this is the case when n — b. 

Theorem 1.1. ^5x5 real matrix is panstochastic if and only if it is a convex combination 
of panstochastic permutation matrices. 
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Because no permutation matrix is a convex combination of other permutation matrices, 
the nxn doubly-stochastic matrices form a generahzed polyhedron whose vertices are the 
n X n permutation matrices [9]. Theorem 1.1 has the following geometric interpretation. 

Corollary 1.1. The set of all real panstochastic 5x5 matrices forms a convex polyhedron 
whose vertices are panstochastic permutation matrices. 

The analogue of Birkhoff's theorem for panstochastic matrices clearly holds if n = 1. 
Our second main result shows that Theorem 1.1 does not generalize for n > 1, n ^ 5. 

Theorem 1.2. If n> 1 and n ^ 5, then there is some nxn panstochastic matrix that is 
not a convex combination of panstochastic permutation matrices. 

We adopt the following notation and terminology: An nxn matrix A over a field is 
panmagic if the sums of the entries along all of the rows, columns, downward diagonals, and 
upward diagonals, both broken and unbroken, of A are equal, in which case the common 
value of the sums is called the magic number of A, denoted IJ,{A). The set of all n x n 
panmagic matrices over F, denoted Pan(n, F), is a subspace of gl(n, F), the vector space 
of all n X n matrices over F. We index the rows and columns of an n x n matrix by 
the elements of ri„ = {0, 1, ... ,n — 1}. We index the diagonals in the following way: for 
k & Qn, the kth upward diagonal contains the {i,j) entry whenever i + j = k (mod n), 
and the A;th downward diagonal contains the {i,j) entry whenever i — j = k (mod n). 

2. Panmagic Permutations. 

Let Sym(ri^) be the group of all permutations on Qn, and let F be a field. For tt G 
Sym(f2„), let P-,^ denote the corresponding permutation matrix: the {i,j) entry of F^- is 
1 if i = 7r(j) and otherwise. We say tt e Sym(O^) is a panmagic permutation of degree 
n if Fn- e Pan(n, F). Denote by the set of all panmagic permutations of degree n. 
These definitions are independent of F; see Lemma 2.1. In particular, if F = R, then 
TT e Sym(f2„) is panmagic if and only if F^. is panstochastic. 

If TT G Sym(On), then the sum of the entries along any row or column of Pt^ is equal to 
1. The sums of the entries along the diagonals of F^. are all equal to 1 if and only if the 
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congruences 7r(j) — j = k (mod n) and Tv{j) + j = k (mod n) have unique solutions for 
j e r^n whenever /c e fi^. We record this observation as 

Lemma 2.1. Let tt e Sym(f2„). T/ien tt is panmagic if and only if there exist X, p & 
Sym(ri„) such that 

7r(j) - j = A(j) (mod n) and 7r(j) +j = p{j) (mod n) 

/or all j E Qn- 

Let a be an integer that is relatively prime to n, and let b be an arbitrary integer. Denote 
by TTax+b the permutation of that sends j to aj+b (mod n). Such a permutation is called 
an affine permutation of 0„. If tt = Tiax+b is an affine permutation, then we sometimes 
denote P,r by Pax+b- Let A„ be the set of all affine panmagic permutations of degree n. 
The following is an immediate consequence of Lemma 2.1. 

Lemma 2.2. Suppose a is relatively prime to n. Then the affine permutation iTax+b is 
panmagic if and only if a + 1 and a — 1 are also relatively prime to n. 

Squaring the congruences in Lemma 2.1 and summing over j gives the following result; 
see [3], [6], and [7]. 

Theorem 2.1. There is some panmagic permutation of degree n if and only ifgcd{n, 6) = 
1. 

This result is already enough to prove one case of Theorem 1.2. If gcd(n, 6) > 1, then 
the panstochastic matrix with all entries equal to 1/n is not a convex combination of 
panmagic permutation matrices. We can therefore rule out this case. 

We next determine which n satisfy n„ = A„. Let n = 5. If tt e lis and 7r(0) = b, 
then the bth entry of the zeroth column of must be 1. This forces all other entries in 
the zeroth column, the bth row, the bth upward diagonal, and the (n — b)th downward 
diagonal to be 0. Consideration of the remaining eight entries shows that either n = TT2x+b 
or TT = 7r3a;+6- Therefore every panmagic permutation of degree 5 is affine: II5 = A5. 
Similar reasoning shows that II7 = Ay and IIn = An. On the other hand, 

Iln 7^ A^ whenever n > 13 and gcd(n, 6) = 1. 
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Indeed, Bruen and Dixon have constructed non-afl&ne panmagic permutations of degree n 
whenever n is prime and n > 13 [2]. See also [7], in which non-affine panmagic permutations 
are constructed whenever gcd(n, 6) = 1, n is composite, and n is not square- free. We 
construct a variation of Bruen and Dixon's example that applies when n is composite and 
gcd(ri, 6) = 1, as follows. Let p be a prime divisor of n. Define S = {k e fin \ k = 
(mod p) } and T = \ S. The permutation n : defined by 

( 2x if X e S, 
^ ^ 1 3a; if a; e T, 

is then panmagic and non-affine. 

3. Kronecker and Wreath Products. 

In this section we consider methods for constructing new panmagic matrices from ex- 
isting ones. Assume that m and n are positive integers and F is a field. Suppose 
^0, • • • 5 An-i e gl(m, F) and B e gl(n, F). Define a matrix {Aq, . . . , An-i)lB e gl(mn, F) 
as follows: if Ag = Taf .1 . . and B = \br s]^ .^n i then the (in + r,jn + s) entry of 
{Ao,...,An-i) IB is 

{{Aq, . . . , An-l) I B).^^^^.^^^ = a- J hr,s 

for z, j e Q,mi r, s e f]n- Notice that if ^ e gl("^, -F), then {A, . . . , A) I B = Ak B, where 
Ak B is the Kronecker product of A and B [4, p. 250] . The next result can be deduced 
from [10, Theorem 2.7 ]. 

Theorem 3.1. Suppose Ao,...A^-i,G Pan(m, F) satisfy ij.{Aq) = ■■■ = ^(An-i), and 
let B e Pan(n, F). Then {Aq, An-i) I B e Pan(mn, F) and {{Aq, An-i) I B) = 
^i{Ao)^i{B). 

Corollary 3.1. If A E Pan(m, F) and B e Pan(n, F), then AxB E Pan(mn, F) and 
li{A xB) = ii{A)ii{B). 

Observe that if Ao,...,A„_i G Sym(OTO,) and p G Sym(On), then the matrix = 
(Fxg, . . . , Fx^_ J ^ Fp is a permutation matrix: the underlying permutation n G Sym(n^„) 
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is given by 

7r(jn + s) = Xs{j)n + p{s) for j e Qrn,s e fin- 
Denote this permutation tt by (Aq, . . . , A„_i) I p. The collection of all such permutations is 
a subgroup of Sym(ri^„) that is isomorphic to the wreath product Sym(riTO) I Sym(ri^). 

Theorem 3.2. // Aq, . . . , ^n-i G Sym(fi^) and p G Sym(n„), then (Aq, . . . , A„_i) I p G 
Hmn if and only if Xq, . . . , A^-i G and p G H^. 

Proof Put TT = (Ao, . . . , An-i) I p. If Ao,...,An-i G and p G n„, then = 
(Pxq, . . . , ? Pp G Pan(mn, F) by Theorem 3.1, and thus n G Hmn- Conversely, 

suppose p ^ n„. By Lemma 2.1 there is some 5 G {—1,1} such that the mapping 
s I— > p(s) + Ss (mod n) is not surjective from to Qn- Hence the mapping x i— > 
%{x) + dx (mod mn) is not surjective from to because 7r(jn + s) + d{jn + s) = 
p{s) + Ss (mod n), and so tt ^ ^mn- On the other hand, suppose A^ ^ for some 
s G fin, so j I— > As(j) + 5j (mod m) is not injective for some 5 G {—1, 1}. For this s and 5, 
the mapping j i-^ 7r(jn + s) + 5(jn + s) (mod mn) is not injective from to rirnn; cind 
hence the mapping x i— > 7r(a;) (mod mn) is not injective from Qmn to Omn- Therefore 

71" ^ n^„. □ 

For A G Sym(nTO), p G Sym(n„), define A x p = (A, . . . , A) ? p G Sym(n^„). Thus 

(A X p)(jn + s) = A(j> + p(s) 

for j G r^rn, s G r^n- The following result is a special case of Theorem 3.2. 

Corollary 3.2. If X E Sy'm{Qjn) and p G Sym(n„) are panmagic, then Xx p is panmagic. 

For A = [ojj] ■ ^.^^2^ G gl(n, F), define supp^ = { (ij) G x ^^n | 7^ }. The next 
result is used in the proof of Lemma 4.2. 

Theorem 3.3. Suppose A G gl(m, F), p G Sym(f2„), tt G Sym.{Qmn)> and suppP^r Q 
supp(AxPp). Then tt = (Aq, . . . , A„-i) Ip for some Ao,...,An-i G Sym(O^), and 
suppPxs Q supped for all s G iln- 



Proof. Let A = [aij]- . Suppose i,j e Qm, r, s G 0„ and Tv{jn + s) = in + r. Thus the 
{in + r, jn + s) entry of P^- is equal to 1, so the {in + r, jn + s) entry oi Ax Pp is nonzero. 
Hence aij{Pp)r,s 7^ 0, and so p{s) = r. Thus there are unique mappings Ao,...,An-i 
from to Qm such that 7r{jn + s) = Xs{j)n + p{s). Since tt is injective, each Xg must be 
injective, and hence Aq, . . . , A^_i G Sym(O^). 

Fix s e flrnj and put A = A^. If the {i,j) entry of Px is nonzero, then A(j) = i, 
so 7r{jn + s) = Xs{j)n + p{s) = in + p{s)- Thus the {in + p{s),jn + s) entry of P-^ is 
1, so the {in + p{s),jn + s) entry of ^ x Pp is nonzero, and hence aij ^ 0. Therefore 
supp Px C supp A. □ 

4. Proofs of Main Theorems. 

Proof of Theorem 1.1. Clearly any convex combination of panmagic permutation matrices 
is panstochastic. Conversely, suppose A G gl(5, M) is panstochastic. According to the 
discussion immediately following Theorem 2.1, lis = A5 = {■n2x+c | c G } U {iTzx+d \ 
d & Q,^}. In [12], Thompson showed that the matrices {P^- | tt G lis} span Pan(5, R). 

Thus there are real numbers etc and (3d such that 

4 4 

c=0 d=0 

'ao + Po "3 + 132 oil + P4 "4 + f3i 012 + /33 
ai + /9i 014 + /^s q;2 + /3o "0 + as + 13a 

0:2 + /52 CKo + /54 0:3 + /?i 0:1+ /?3 04 + /?o 
as + /5s CKl + I3q 0:4 + /?2 0:2 + /54 CKO + /^l 

.a4 + P4 a2 + Pi ao + Ps a^ + Po OL1+P2- 
If a cyclic permutation of rows or columns is applied to a matrix of the form P2a;+c {P^x+d-, 
respectively), the result is another matrix of the form P2x+c {Pzx-\-di respectively). Hence 
we can assume without loss of generality that uq — minjctc} and /3o = min{/3d}. Also, 
ao + /^o ^ and Ylic^c + Ylidl^d = 1 because A is panstochastic. If both ckq and Pq 
are nonnegative, then (1) expresses ^ as a convex combination of panmagic permutation 
matrices. If /3o < 0, then 

4 4 

A = ^{ac + Po)P2x+c + J2^Pd- Po)Pzx+d 

c=0 d=0 



(1) 



- 7 - 



is a representation of the required form, and a similar expression can be obtained if ckq < 
0. □ 

We have aheady dispensed with the case gcd(n, 6) > 1 in Theorem 1.2. The remainder 
of the proof uses the following two lemmas. 

Lemma 4.1. If n > 1 and gcd(n, 30) = 1, then there is some panstochastic n x n matrix 
that is not a convex combination of panmagic permutation matrices. 

Proof. The 7x7 matrix 

■0 1 1 0- 

1 1 

1 1 

2 

1 1 

1 1 

.0 1 1 OJ 

is panstochastic. Suppose tt = iTax+b G Ay = Ily and supp^Tr ^ suppA. Then b = 7r(0) G 
{1,2} and 7r(3) — 3a + b = 3. If 6 = 1, then a = 3, so 7r(2) = 0, which contradicts 
suppPtt Q supp^. On the other hand, if 6 = 2, then a = 5, so 7r(l) = 0, which again gives 
a contradiction. Thus the assertion of the lemma holds when n = 7. 

For the rest of the proof we suppose n > 7 and gcd(n, 30) = 1, so n > 11. Observe 
that n2x+i and tt2x-4 are affine panmagic permutations of degree n by Lemma 2.2. Define 
^0 = (l/2)(-P2a;+i + P2x-4)i SO is panstochastic and takes the form 



2 



Ao = - 
" 2 



10* 

1 * 
1 * 
10 0* 



Consider the n x n matrix 



1 
2 




-1 
1 





1 



-1 





-1 




1 







1 
-1 
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with zero entries except in rows through 3 and columns through 3. Observe that 
Ai e Pan(n, R) and iJ,{Ai) — 0. Define A — Aq + Ai, so A is panstochastic and has the 
form 



10 0* 

1 * 

1 * 
10* 



Assume that A is a convex combination of panmagic permutation matrices. Then there 
exists some tt e n„ such that 



7r(0) — 2 and suppP-^- C suppA. 

We show that the existence of such a permutation tt leads to a contradiction. 
Notice that from (2) we have 



(2) 



7r(l) e {0, n - 2} , 7r(2) e {3, 5} , 7r(3) G {1, 7} . 

Put Un = {0, 1,2,3}. If £ e n„, then either 7r(£) = 21+1 (mod n) or 7r(£) 

2£-4 (mod n) by (2). Define subsets J, K ofQ n as follows: 

{j eTln \ Tv{j) = 2j + l (mod n) }, 

K= {keUn I 7^{k) = 2k-4: (mod n) }. 
Since tt is panmagic, we have 

7r(i) — i ^ 7r(j) — j (mod n) whenever i, j e and i ^ j 



(3) 



(4) 



by Lemma 2.1. Moreover, 

if J e J and j < n — 5, then j + 5 e J. 

Indeed, if j e J and j + 5 e K, then 

7r(j)-j = (2j + l)-j= (2(i + 5)-4)-(i + 5) 

= "^(j + 5) — {j + 5) (mod n). 



(5) 



contradicting (4). It follows that 

iik e K andk>9, then k-5eK. (6) 

since J and K form a partition of Qn- 

The following statements (7) - (10) all follow from (4). 

If n - 1 e J, then 4 e J. (7) 

If n - 4 e J, then n-2e K. (8) 

If n - 3 e J, then 8 e J. (9) 

If 7 e then 5 e J. (10) 

For (7), suppose n — 1 E J and 4 E K. Then 7r(n — 1) — (n — 1) = = 7r(4) — 4 (mod n), a 
contradiction. For (8), suppose n — 2 e J and n — 4 e J. From the definition of J we have 
7r(n — 2) — (n — 2) = n — 1 (mod n) and 7r(n — 4) — (n — 4) = n — 3 (mod n). However, 
by (3), 7r(l) — 1 = n — 1 (mod n) or 7r(l) — 1 = n — 3 (mod n), which is a contradiction. 
The proofs of (9) and (10) are similar to that of (8), using the values of 7r(3) and 7r(2), 
respectively. 
We have 

j E J whenever 4 < j < n and j = 1 (mod 5). (11) 

Indeed, by (5), it is enough to show 6 e J. For this, observe 7r(0) — = 2, so 7r(6) — 6 7^ 2 
by (4), and hence 6 ^ K. 
In addition, we have 

^ e K. (12) 
Indeed, if (n + l)/2 e J, then tt ((n + l)/2) = n + 2 = 2 = 7r(0) (mod n), a contradiction. 

We finish the argument by cases, according to the congruence class of n modulo 5. 
Case 1: n = 1 (mod 5). In this case (n + l)/2 = 1 (mod 5), and thus (n + l)/2 e J by 
(11). This contradicts (12). 

Case 2: n = 2 (mod 5). In this case, n — 1 = 1 (mod 5), and thus n — 1 e J by (11). 

By (7), 4 G J, and by (5), j E J whenever 4 < j < n and j = 4 (mod 5). Since 
(n + l)/2 = 4 (mod 5), (n + l)/2 e J, contradicting (12). 
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Case 3: n = 3 (mod 5). Since {n+l)/2 = 2 (mod 5), it follows from (12) and (6) that 

7 & K, and thus 7r(7) = 10 (mod n). However, if j = (n + 9)/2, then 4 < j < n and 
j = I (mod 5), so j E J by (11), and hence 7r{j) = 2j + 1 = 10 (mod n), which is a 
contradiction. 

Case 4-' n = 4 (mod 5). In this case, n — 3 = 1 (mod 5), and so n — 3 G J by (11). Hence 

8 e J by (9), and hence n — 1 e J by (5) since n — 1 = 8 (mod 5). Thus 4 e J by (7), 
and so j E J whenever 4 < j < n and j = 4 (mod 5) by (5). Since n — 4 = 4 (mod 5), 
we conclude n — 4 e J. By (8), n — 2 e and since n — 2 = 2 (mod 5), we have 7 e K 
by (6). By (10), 5 G J, and by (5), j E J whenever 4 < j < n and j = (mod 5). Since 
(n + l)/2 = (mod 5), we have (n + l)/2 G J, contradicting (12). 

We have arrived at a contradiction in each case. □ 

Lemma 4.2. Suppose there is some mxm panstochastic matrix that is not a convex combi- 
nation of panmagic permutation matrices and that n is a positive integer with gcd(n, 6) = 1. 
Then there is some (mn) x (mn) panstochastic matrix that is not a convex combination of 
panmagic permutation matrices. 

Proof. There exists a panmagic permutation p of degree n, say p = tt2x- Among all 
panstochastic mxm matrices that are not convex combinations of panmagic permutation 
matrices, choose one, say A, with the maximum number of zero entries. The matrix 
B — Ax Pp is then a panstochastic (mn) x (mn) matrix by Corollary 3.1. Assume that B 
is a convex combination of panmagic permutation matrices. Then there is some panmagic 
permutation tt of degree mn such that suppPvr C suppS. By Theorem 3.3, there are 
Ao, . . . A^_i G Sym(O^) such that tt = (Aq, . . . , A„_i) I p and suppP^^ ^ supp A for all 
s G Then Aq, . . . , A„_i are panmagic by Theorem 3.2. Fix some s G i^m and put 

A = As, so suppP\ C supp A. Define 

a = mill {ttij I A(j) = i}. 

We cannot have a — 1, for otherwise A = Px is a. panmagic permutation matrix. Therefore 
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< a < 1, and so the matrix 

1 — a 1 — a 



is a panstochastic mxm matrix with a greater number of zero entries than A. By the choice 
of A, C is a, convex combination of panmagic permutation matrices, and it follows that 
A = aP\ + (1 — a)C is also a convex combination of panmagic permutation matrices, so we 
have a contradiction. Therefore B is not a convex combination of panmagic permutation 
matrices. □ 



Proof of Theorem 1.2. By Lemma 4.1 and Lemma 4.2, it suffices to prove that there is some 
25 X 25 panstochastic matrix that is not a convex combination of panmagic permutation 
matrices. The matrix 
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1. 



is panstochastic; it can be obtained by averaging the matrices for t^2x+i and a non-affine 
panmagic permutation and then adjusting the entries in rows through 3 and columns 
through 3, as in the proof of Lemma 4.1. Suppose ^ is a convex combination of panmagic 
permutation matrices. Thus there is some panmagic permutation tt on such that 
7r(0) = 2 and suppP^r Q supped. Since 7r(0) = 2, we cannot have 7r(13) = 2, and hence 
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7r(13) = 18. Thus n{21) 7^ 18, so 7r(21) = 4. Also, 7r(0) ^ 17, so 7r(8) = 17. Therefore 
7r(8)+8 = = 7r(21)+21 (mod 25), so n is not panmagic by Lemma 2.1, and a contradiction 
is reached. □ 
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